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Abstract 

- - - ,We extend to large lattices the work of a previous investigation of the phase diagram of the anisotropic five-dimensional 
I SU(2) Yang-Mills model using Monte Carlo simulations in the regime where the lattice spacing in the fifth dimension is 
Jarger than in the other four dimensions. We find a first order phase transition between the confining and deconfining 
fSl phase at the anisotropic parameter point /34 = 2.60 which was previously claimed to be the critical point at which the 
j^prder of the transition changes from first to second. We conclude that large lattices are required to establish the first 
d order nature of this line of transitions and consequently that the scenario of dimensional reduction of the five-dimensional 
^1 Itheory to a continuum four-dimensional theory via the existence of the so-called "layer phase" is unpromising. 
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1. Introduction 
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Since the idea of introducing extra dimensions to give 
a better understanding of the hierarchy between weak and 
Planck scales was introduced more than a decade ago in [ij , 
many models have been developed to show how extra di- 
mensional theories might solve these problems and yet ex- 
liibit a connection to the four-dimensional world that we 
observe. Well-known models that lead to dimensional re- 
'duction are those of Randall-Sundrum (RS) ADD P, 

0,0], Dvali-Shifman mechanism (DS ) B a nd compactifi- 
'cation of the extra dimensions 0, lol llOf and since the 
.work of 11| they have gained phenomenological interest. 

The first three models envisage the four-dimensional 
.world as layers that exist in the extra dimensions and are 
'decoupled from each other. All the particles that we cur- 
.rently observe can only propagate in one layer, with the 
'exception of the graviton that can interact between the 
layers. The RS model deals with warped space-time, while 
the ADD model (D-brane models) and DS mechanism con- 
sider the extra dimensions as fiat. The compactification of 
the extra dimension is based on Kaluza and Klein's at- 
tempt to unify fundamental forces in the early 20th cen- 
tury [ll[i3. 

Since higher-dimensional non-abelian models are non- 
renormalizable, it is necessary to define an effective theory 
by defining parameter values at which there is a separa- 
tion of scales between the cutoff A and low energy physics. 
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The first study of the localization of fields on branes using 
non-perturbative methods came from Fu and Nielsen lj| . 
They propose that, when there is a Z?-dimensional lat- 
tice, a d-dimcnsional layer phase [D — n + d) can be 
formed if the nearest-neighbour gauge couplings of the 
d-dimensional sublattice are different from the other n 
dimensions. Then, particles and gauge fields can travel 
within the d-dimensional layer phase, but they exhibit a 
kind of confinement when they try to propagate in any of 
the n extra dimensions. The existence of the layer phase 
was shown in the five-dimensional U(l) gauge theory using 
Monte Carlo simulations in [l5, 16, 17 1. 

When the lattice SU(2) Yang-Mills model is considered 
in five dimensions, the so called "layer" phase is believed to 
exist when the anisotropy in the lattice couplings is such 
that the lattice spacing in the extra dimension is larger 
than the one in the usual four dimensions, i.e. > a^. 
In this layer phase one expects a zero string tension in the 
usual four-dimensional directions, but a non-zero string 
tension along the extra dimension. 

The phase diagram of the five-dimensional non-abelian 
gauge theory was first studied in 1979 by Creutz [2l|, who 
showed a phase transition between the confined and the 
deconfined phase of the model with only isotropic cou- 
plings. During the last decade, the dimensional reduction 
from the five-dimensional theory to four dimensions, by 
implementing the model on a torus, was investigated on 
an anisotropic lattice using Monte Carlo techniques 0, H, 
These papers show that, when the model is di- 
mensionally reduced via compactification, the phase tran- 
sition changes its nature to second order and belongs to 
the same universality class as the four-dimensional Ising 
Model. In 0, i, [I^ simulations were done in the region 
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where 04 > 05 , whereas in [9| the simulations were carried 
out in the anisotropy region which is of interest in the case 
of the layer phase, i.e. 04 < 05. 

The existence of the layer phase was investigated by 
mean-field approximation on an anisotropic lattice with 
periodic boundary conditions 18, [l^ and it was shown 
that the planes transverse to the extra dimension were de- 
coupled from each other. An investigation of the existence 
of this layer phase using Monte Carlo techniques was at- 
tempted by Farakos et al. [2^ . They claim that the transi- 
tion between the five-dimensional Coulombic (deconfined) 
phase and the strong-coupling (confined) phase changes 
its order from first to second, implying that in the fifth 
dimension the layer phase exists. This opens a possibility 
of defining a continuum four-dimensional field theory. 

The motivation for studying this model further came 
mostly from results in Q, which indicate that the lattice 
volumes used in the previous Monte Carlo simulations (20j , 
were too small to show the correct order of the phase tran- 
sition. In this letter we extend this investigation to larger 
volumes. 

This letter is structured as following: In Section 2 we 
set up the lattice model and the observablcs that were 
measured. In Section 3 we give details for our simulations 
and we present our results, and in Section 4 we give a brief 
conclusion. 



2. The Model 

2.1. Anisotropic Action 

We investigate the anisotropic SU(2) Yang-Mills gauge 
theory in five dimensions, whose action in the continuum 
is given by 



Se = I d^x I dx^^ T'rF^jpf 
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where M,N ~ 1...5 and Fmn = ^a/Aat — djqAu + 
i{AM,AN\ with Am = g^AI/T". 

On the lattice the action of the model becomes 



X l<p<i/<4 



+^5^ E (1-^ TrC/^5(x-)) y.,.=l. 
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where \J (x) represents the oriented plaquette along space- 
time directions given by 

Uf,^{x) = Uf,{x)U^{x + flaiplix + i>a4)Ul{x) (3) 

and U^^lx) represents the plaquette formed when one of 
the directions is the extra-dimensional one, given by 

= U^{x)Uz{x + (la^Wlix + ^^a^)Ul{x). (4) 

where = exp(if/5a4A^) and U5 = eyL-piigc^a^Ac,) are the 
gauge links and 04 is the lattice spacing in the temporal 



and three spatial directions and 05 is the lattice spacing 
in the extra direction. 

The anisotropy parameter on the lattice is character- 
ized by 7 which is given by 



7 = 

and at classical level this is given by 

a4 
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2.2. Observablcs 

In order to investigate the phase diagram of the model 
we use the following observablcs: 



• Average Plaquette in the extra dimension a;5, P5 

■J2J2^TiU,,ix))) (7) 



and its susceptibility 



(8) 



where V is lattice volume given hyV^ Lt x L| x L5, 
with Lt, Ls and L5 the size of the temporal, spatial 
and extra dimension respectively. 

Temporal Polyakov Loop. This can be measured on 
the whole lattice given by 



Poly, 



(LT-l)a4 

^ Tr n t^i(^) 

x,X5 a;i— 



(9) 



Since, in the layer phase, each layer is uncorrelated, 
the Polyakov loop may be measured in one layer and 
so we also compute 



Poly^(a;5) 
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{LT-l)a4, 

5^ Tr n U,{x)\,, 

X xi=0 



(10) 



We define the Polyakov loop susceptibilities as 
V 



XPoly, 



(Poly^ - (Poly^)' 



(11) 



XPoly (2^5) = i|((Poly^(x5)2 - (Poly^(x5))2) 

(12) 

The expected behaviour of the plaquette for a first order 
phase transition is to show hysteresis in the expectation 
value and a divergence in the susceptibility at the critical 
point. The temporal Polyakov loop is expected to have a 
zero expectation value in the strong phase, i.e. to fluctu- 
ate around zero and a non-zero expectation value in the 
deconfining phase, i.e. to show a two-peak structure. 
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3. Results from Lattice Simulations 

Our model was implemented on the lattice, using the 
Kennedy-Pendleton Heat-Bath algorithm combined 
with overrelaxation updates Q ■ Specifically, we took one 
heat-bath measurement every Ls/2 overrelaxation steps. 
The autocorrelation that arises was taken into account in 
our analysis. The number of measurements varied between 
100,000 and 200,000 at each set of points in our param- 
eter space (/34, /Js) that were investigated. Measurements 
were taken starting from either random SU(2) matrices 
(hot configurations) or by setting all the SU(2) matrices 
to the identity matrix (cold configuration). The lattice 
volumes were 16^, 20"^ x 8 and 24^ x 8. The former was 
the largest investigated by Farakos and Vrentzos [2^ and 
was included here as a check with their results. The bigger 
volumes were used, since Q showed that there is a mini- 
mum size of the spatial/temporal and the extra direction 
in order to see a clear first order phase transition. Since, 
16^ is below this minimum size, we simulated bigger vol- 
umes to investigate the order of the phase transition. We 
reduced the size of the extra dimension to L5 = 8 to save 
compute time, since the lattice spacing in the extra di- 
mension is much larger for high /S^ values than the lattice 
spacing in the other directions (as shown in Q), so the 
system remains five dimensional. 
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Figure 1: A sketch of the phase diagram of the anisotropic SU(2) 
Yang-Mills model. The dashed blue line denotes the isotropic case 
_7 = 1. The region above this line was previously investigated in 0, 
[a, |10|| and the region below in [91. 12011 . The dashed-dotted green line 
appears when the extra dimension is compactified 0, 0, [Toll . When 
no compactification is involved, there is a bulk phase transition which 
is shown in the figure as a red solid line. It was shown to exist up 
to /J4 = 2.50 in [31. In this work we extend the range of this line 
up to = 2.60 with no evidence that this line will not continue for 
larger values of /34. For ^4 > 2.60 the idea of the existence of the 
layer phase arises. 

First, we did a scan in the parameter space (/34, jS^) 
using small lattices to identify the first order phase tran- 
sition that was shown in previous work up to (3^ ~ 2.50. 
The phase diagram of the model is shown schematically in 
Fig. [1] The layer phase was previously claimed to exist at 



large /S^ and small P^, as shown. Our point of interest is 
(34 = 2.60 on the line of transition, which was claimed to 
be the critical point at which the transition changes from 
first to second order in The critical point in /?5 was 
found by implementing the model on a lattice of volume 
16^. At this volume we were able to do a wide scan by 
investigating a sufficient number of different (3^ to identify 
the critical point. Even though it does not show any clear 
evidence of first order phase transition in terms of a two- 
state signal, by looking at the susceptibility it looks like it 
has a divergence at the critical point (Fig. The criti- 
cal (3^ point was found to be f3^ = 0.8437(5) which agrees 
within error with the value found in (20j . 
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Figure 2: The susceptibility of the plaquette in the extra dimension, 
P5 for V = 16^^ keeping ^4 fixed at 2.60 and varying /Ss . The critical 
point is the point at which the susceptibility gains its maximum 
value. 

For the investigation of the phase transition on the 
larger lattices we focused on the critical region which was 
estimated to be between /^s = 0.843 and ^5 = 0.8445, 
based on the critical value found for the 16^ lattice. As can 
be seen from Fig. [31 the plaquette moves to the right as we 
go to higher values of /3^ and the temporal Polyakov loop 
is zero for the point /3^ = 0.843, which is in the confining 
phase and has a two-peak structure for the point /Ss = 
0.8445, which is in the deconfining phase, as expected. We 
also checked how the temporal Polyakov loop is distributed 
when all the xs-slices are considered independently (Eq. 
[TU)) and we could see that each fiuctuates around zero. 
Also, we confirmed that the critical point was included 
in this region by observing that, for one point that lies 
in between these values, either a clear two-state signal or 
large fiuctuations between two values in the average value 
of the plaquette were present, as can be seen in Figures 2] 
and[5l 

The points that were investigated are (3^ = 0.843, 0.8435, 
0.844, 0.8445. For the lattice volume of 20^ x 8, we cannot 
distinguish the two states, since by starting either from hot 
or cold configurations, the expectation of the plaquette is 
the same. However, as shown in Fig. 01 the distribution is 
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Figure 3: The histograms for the plaquette in the extra dimension 
for V = 20* X 8, fixed jS^ = 2.60 and for two different values of 
/Ss: jSs = 0.843(top) and ^5 = 0.8445 (bottom) are shown on the 
left. We can see that the peak moves towards the right as we go to 
higher values of /^s. The corresponding histograms for the temporal 
Polaykov loop are shown on the right. We can see that for /Js = 0.843 
it has a zero expectation value whereas for fis = 0.8445 it shows a 
two-peak structure. 



not Gaussian. This is the first hint of the existence of a 
first order phase transition. The larger volume of 24* x 8 
shows a clear two-state signal as can be seen from the his- 
togram in Fig. [5] The overlap that appears between the 
histograms starting from hot or cold configurations is due 
to the fact that the two vacua of the potential energy are 
not so deep and thus the system fluctuates between them. 
For a check, we implemented one single point (/Ss = 0.844) 
on a 32'' x 8 lattice and we can see in Fig. \6\ that the two 
states are now separated by a wide gap, and there is no 
evidence of tunnelling from one to the other, i.e. it stays in 
the phase in which it first equilibrated, depending on the 
initial configuration. This is also an indication that the 
extrapolation to the thermodynamic limit must be based 
on sufficiently large lattices. 

We note that the critical point was not estimated pre- 
cisely, because reweighting techniques were not trustwor- 
thy for the large volumes and ensemble sizes that were 
used in this work due to limited statistics. 

The code was written using QDP++ [l^l and run on 
GPUs using QDP-JIT [25j. An estimate of the compute 



time required for the investigation of a single point on 
an NVIDIA Tcsla C2070 Computing Processor (GPU) for 
the volumes used in this work for a set of 100,000 mea- 
surements with Ls/2 ovcrrclaxation steps and a heatbath 
update each time is shown in Table[T] The GPUs were pro- 
vided by the Particle Physics Theory Group at the Univer- 
sity of Edinburgh and the Edinburgh Compute and Data 
Facility. The single point oiV^ 32** x 8 would have taken 
two months on GPUs and so was simulated using STFCs 
DiRAC facilities in Edinburgh. 
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Figure 4: Histograms of the average plaquette in the extra dimension, 
P5 starting from both cold and hot configurations for V = 20* X 8, 
04 = 2.60 and /35 = 0.8435. We can see that since this point is very 
close to the critical one, the plaquette fiuctuates between the two 
vacua and thus the distribution is not Gaussian anymore. 
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Figure 5: Histograms of the average plaquette in the extra dimension, 
P5 starting from both cold and hot configurations for V = 24* X 8, 
j34 = 2.60 and f}^ = 0.8435. Here, we can see that the distributions, 
starting from either cold or hot, build up as two Gaussian distribu- 
tions, one for each vacuum that the system equilibrates to. 
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Figure 6: Histograms of the average plaquette in the extra dimension, 
P5 starting from both cold and hot configurations for V = 32* X 
8, /34 = 2.60 and fis = 0.844. It is clear that a first order phase 
transition is present since starting from different configurations, the 
system equilibrates in different states with no tunnelling between 
them. 
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Lattice Volume 


Compute time (hours) 


16 X 16 X 16 X 16 X 16 
20 X 20 X 20 X 20 X 8 
24 X 24 X 24 X 24 X 8 


190 
250 
620 



Table 1: Estimated compute time required on an NVIDIA Tesla 
C2070 Computing Processor for 100,000 measurements for a single 
point in the parameter space (/34,/35). 



4. Conclusions 

In this work, we extended the Monte Carlo investiga- 
tion of the phase diagram of the anisotropic SU(2) Yang 
Mills model in five dimensions when the lattice spacing in 
the extra dimension is larger than that in the four other 
dimensions (7 < 1). We showed that, up to /34 = 2.60, 
there is no evidence of a second order phase transition, 
whereas a clear two-state signal in the average plaquette 
favours a first order phase transition. Based on this result, 
we can claim that the bulk first order phase transition be- 
tween the confining and the deconfining phase continues at 
least up to /34 = 2.60 with no trace of a layer phase. This 
implies that up to this point the continuum limit cannot 
be taken and thus the possibility of a dimensionally re- 
duced five-dimensional effective field theory remains open. 
Even though, based on the work of this letter, we cannot 
exclude a second order transition at higher /34, nothing in 
our study suggests that continuing this investigation on 
even bigger lattices would be worthwhile. 
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